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We investigate analytically and numerically a random-matrix model for m fermions occupying 
li single-particle states with positive parity and £2 single-particle states with negative parity and 
interacting through random two-body forces that conserve parity. The single-particle states are 
completely degenerate and carry no further quantum numbers. We compare spectra of many- 
body states with positive and with negative parity. We show that in the dilute limit defined by 
QQ I m,£i^2 —> oo,m/£i^2 — > 0, ground states with positive and with negative parity occur with equal 

, probability. Differences in the ground-state probabilities are, thus, a finite-size effect and are mainly 

' due to different dimensions of the Hilbert spaces of either parity. 
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PACS numbers: 21.60.Cs,24.60.Lz,21.10.Hw,24.60.Ky 



I. MOTIVATION AND AIM 



■ Johnson et al. [l| observed that in the two-body random ensemble (TBRE) of the nuclear shell model, ground states 
with spin zero occur much more frequently than corresponds to their statistical weight. That observation caused 
considerable theoretical activity (see the reviews [2] and [3]). A similar preponderance for states with positive parity 
was found in Ref. We wish to explore the reason for that preponderance. We focus attention on parity (rather 
than spin) because that quantum number is analytically more easily accessible. We use a model with spinless fermions 

^ ' that interact via random two-body forces. The degenerate single-particle states carry no orbital angular momentum 
^ [ quantum number but have either positive or negative parity. The model is a modified version of EG0E(2), the 
embedded two-body ensemble of Gaussian random matrices Q. We investigate the model by using both, an analytical 
approach and numerical simulations. The analytical approach evaluates traces of powers of the Hamiltonian up to very 
high order and uses results of Refs. and to estimate the position of the ground state. The numerical simulations 
' involve diagonalization of matrices drawn at random from the ensemble and can be done only for Hamiltonian matrices 
f — ■ , of sufficiently small dimension whereas the analytical approach is suited also for large-dimensional matrices. 
On ' To motivate our focus on traces of the Hamiltonian, we recall in Section |TT] how the ground-state energy was 
] estimated in Refs. @ and That method is used and compared with numerical simulations in Section |Vl Prior to 
00 ■ that, we define our model in Section Hill The first and second moments of the Hamiltonian are calculated for both 
' parities in Section IIVI After presenting our numerical results, we investigate our model in the limit of large matrix 
00 , dimension N in Section IVTl We show that for N ^ 00, both the first and the second rnoments of the Hamiltonian have 

■ the same values for either parity. Combining that fact with the well-known result |5| that the shape of the average 
J> , spectrum is asymptotically {N 00) Gaussian, we conclude that ground states of either parity are equally likely. 

In Section IVHI we show that the strong correlations found asymptotically for the first and second moments extend 
^ to higher (but not to all) moments. We discuss the implications of that result for correlations between the spectral 
fluctuation properties of positive- and negative-parity states and show that the result reinforces our conclusions. We 
conclude with a summary and discussion. 



II. SIMPLE ESTIMATE FOR THE GROUND-STATE ENERGY 



To estimate the ground-state energy, we use with proper modifications the method introduced for states with spin 
in Ref. Q and improved in Ref. Q- Let H denote the Hamiltonian of the system, V± the projectors onto states with 
positive and negative parity, respectively, and i?ground(i) the energies of the lowest state with positive or negative 
parity, respectively. We estimate i?ground(±) by writing 

Aground (±) = nTrace(iJ7'±) - r±a± . (1) 

The symbol nTrace stands for the normalized trace (the actual trace divided by the dimension N± of Hilbert space) , 
and the width a is defined as 



cr^ = nTrace(77^7'±) 



(2) 
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In Ref. [g the analogue of Eq. ([T]) was used without the first term on the right-hand side. That term was added in 
Ref. Q ■ It represents the fluctuations of the centroid of the spectrum. Inclusion of that term improves the agreement 
with numerical simulations: the fluctuations of theparameter r are reduced. Equation ([T]) has a simple interpretation: 
shell-model spectra have nearly Gaussian shape [S'l and are, thus, essentially characterized by the centroid and the 
width. The distance of the lowest state from the centroid of the spectrum is given by a multiple r± of the width. In 
the case of spin, the stochastic fluctuations of r were found to be small, so that r can be considered a constant. In 
Ref. |7j|, an explicit expression for r was obtained by fltting the results of numerical calculations. It reads 

r = VO.gginTV -1-0.36 . (3) 

We actually prefer to determine r± by a fit to numerical data. In Section |V] we compare the result with Eq. ([3]). We 
also compare the numerically determined probability of finding a ground state of given parity with predictions derived 
from Eqs. ^ and ©. 



III. MODEL 



We consider a system of m spinless fermions distributed over a set of degenerate single-particle states. There 
are li states of positive parity and ^2 states of negative parity, with associated creation and destruction operators 
a\^, a^^ {ii = l,2,...,i?i) and a\p, {p = 1,2,..., ^2), respectively. The single-particle states carry no further 
quantum numbers. The many-body states of the system have positive (negative) parity if the number 7712 of fermions 
in negative-parity states is even (odd, respectively). The total numbers and of positive- and negative-parity 
states are 

+ / , ^mi+m2,m^m2,cvcn I I I I i 

^- = X! '5™l+m2,m'5m2,odd( ( ^ j • (4) 

mi,m2 V^a/ 

The Hamiltonian is a sum of two-body interactions that conserve parity, 
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The ranges of the summation indices depend in an obvious way on the creation operators and matrix elements on 
which they appear. The two-body matrix elements obey the symmetry relations 
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An ensemble of Hamiltonians is obtained when we consider the matrix elements in Eq. fS]) as Gaussian-distributed 
random variables. We assume that the V^fi^po- are not correlated with the V^^i.pi^i and likewise for the pairs {Vpv-^pa^ 

Xl!^l';p'a') for j = 1,2 and fc = 1,2, and for the pair xl])-pa, -'^^v' p'ct'- matrix elements have zero mean values. 
For the variances, we define pairs of indices a, (3 by writing a — {fJ-iy} and likewise for /3, and have for i = 1, 2 



X^X s' = v'S^^,Sf,p, . (7) 
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The bar denotes the average over the ensemble, and Sap stands for [S^^j^iSy^i — 5^,yi5^^i), etc. The matrix elements 
X^"^^ do not possess any symmetry properties and obey 

X'i^J.paX'^^}^,.^,^, = v'^5^j_^j_>5^i,^5pp'5aa' ■ (8) 

Without loss of generality, we put i;^ = 1 in the sequel. 



IV. CALCULATION OF nTrace(//) AND OF nTrace(7j2). 

These two traces are needed for the evaluation of Eqs. ([T]) and The only non- vanishing contributions to the two 
traces arise from terms in H and in which leave the number of fermions in every single-particle state unchanged. 
These terms are found by using Wick contractions of the creation and annihilation operators in the expressions for 
H and H^. We indicate the omission of all other terms by an arrow. For H we obtain 

H ^\Y1 ^iiU'^iA'"!!^ + ^ 51 ^i^U"2/.ri2;. + Xflppni^n2p . (9) 

/ip 

Here n.i^ is the number operator for state (i/i) with i = 1,2. 

The diagonal element of ni^ni^ taken between one of the states with mi fermions in positive-parity single-particle 
states and TO2 fermions in negative-parity single-particle states vanishes unless both states and (li') are occupied, 
in which case the matrix element equals unity. There are altogether {^~^2) im ) ^^^ch states. We consider separately 
the normalized traces over the positive-parity and the negative-parity many-body states. We recall that V± are the 
projection operators onto the many-body states with positive and negative parity, respectively. We obtain 



nTrace(i/7'+) = ^^^Fjil^^ ^ 6,m+m2,mSm2,cvcn(^ _2 
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nTrace(i/-p_) = --^ V V^^l] V 5„^^+,n^^,n5v^2,odA , ^ , , 



Both traces are seen to depend on the same three uncorrelatcd random variables, 



(10) 



-1 = E ^^iU ; -2 = E ; -3 = E ^^^PP ■ (11) 

/^p 

As sums of uncorrelated random variables with equal Gaussian distributions, zi, Z2, and Z3 have Gaussian distributions 
with mean values zero and second moments — l)/4], [^'2(^2 ~ l)/4], and i'i^2, respectively. Thus, the distribution 

of the traces in Eq. PH)) is completely known. 

The pattern that emerges in Eq. (fTO|) will be seen to apply quite generally to traces of arbitrary powers of H: the 
traces are sums of products. The first factor in each product depends only on the random variables and is the same 
for both parities. The second factor differs for states of positive and states of negative parity but is independent 0/ the 
random variables. That general pattern will be decisive for our understanding of the preponderance of ground states 
with positive parity. 

We turn to Trace (i?^). The following terms yield non-zero contributions: the square of the first term on the 
right-hand side of Eq. ([5]) , the square of the second term, the product of the first and the second term, the square of 
the third term, and the square of the fourth term. We consider these terms in turn. 
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In the square of the first term, there appear the two matrix elements V^^'' with their associated creation and 
annihilation operators. Wick contraction is possible in three different ways: (i) We contract the two creation and the 
two annihilation operators associated with the same matrix element. That is the same procedure as used in formula ([9|) 
and yields a total of 4 contraction patterns, (ii) We contract one of the two creation operators associated with the 
first matrix element with an annihilation operator associated with the same matrix element, and the other with an 
annihilation operator associated with the second matrix element. That yields a total of 16 contraction patterns, (iii) 
We contract the two creation operators associated with the first matrix element with the two annihilation operators 
associated with the second matrix element. That yields a total of 4 contraction patterns. It is straightforward to 
check that because of the fermionic anticommutation rules and the symmetry properties ([6]), the different contraction 
patterns in each of the three groups yield identical results. For the square of the second term on the right-hand side 
of Eq. (O, these considerations apply likewise. For the product of the first and second term, only the contraction 
patterns used in formula ^ are possible. In the square of the third term on the right-hand side of Eq. ([5]), only the 
product of the two terms in round brackets gives a non-vanishing contribution, with obvious contraction patterns. In 
the square of the fourth term, there are the same three possibilities as in the square of the first term. Altogether this 
yields 
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Before working out the trace of this expression, it is useful to rearrange it in such a way that in all summations no 
two summation indices take the same values. This yields 
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In calculating the trace, we observe that the number of nonequal summation indices in the terms on the right-hand 
side of Eq. (fT5)) determines the weight factors. The result is 

nXrace(-/? ^+) ^ ^ f^mi +m2 ,mf^m2 ,cvcn 
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For nTrace(iJ^7'_) we find exactly the same expression except tliat the second Kronecker delta in the first line on the 
right-hand side of Eq. (fT4|) is replaced by 5^2, odd, and that N+ is replaced everywhere by N^. 

As in the case of TraceiJ, the trace of is a sum of terms each of which is the product of two factors. One 
factor depends only on the random variables and is the same for both parities. The distribution of these factors can 
be worked out. That is not done here. Some of the factors are correlated with each other. The other factor is a 
weight factor that is a sum over products of binomial factors. It does not depend on the random variables and is not 
obviously the same for the two parities. Our result would not apply in the case of states with spin where the linear 
or bilinear forms containing the random variables will depend on the total spin. Assuming that Eqs. ([1]) to ([3]) hold, 
we conclude that a preponderance of ground states with even parity - if it exists - can have only one of two causes: it 
may be due to differences between the non-statistical weight factors, or to differences in the scale factors r+ and r_ . 
(We recall that according to Eq. ([3]), the latter depend on the matrix dimensions N±.) 



V. NUMERICAL RESULTS 



For a test of Eq. ([T]) , we perform numerical simulations. To this purpose, we consider several systems that differ in the 
parameters ii, £2, and m. For each set of parameters, we set up the matrix corresponding to the Hamiltonian ([5|) in a 
space of Slater determinants. The Gaussian-distributed two-body matrix elements are computed by a pseudo-random 
number generator, and the ground-state energies i?g,.ound(i) are obtained from a numerical diagonalization of the 
Hamiltonian matrix. For the largest-dimensional matrices, we employ the Arpack package @ in the diagonalization. 
In addition to the ground-state energy we also compute the normalized traces nTrace(i?'^7'±) for fc = 1,2. Our 
ensemble consists of 100 random Hamiltonians for each set of parameters £1, £2-, and m, and we record the ground- 
state energies £^ground(±) and the first two moments nTrace(i?'^7-'±) (with fc = 1, 2) of the parity-projected Hamiltonian 
for each member of the ensemble. We employ Eq. ^ and determine the scale factors r± that relate the ground-state 
energy to the first and second moment of the Hamiltonian by fit. An example is shown for the set of parameters 
m = f 1 = £2 = 9 in Fig. [TJ The results obtained for the scale factors (with their rms variances) are shown in Table [D 
The table also shows the probability p+ that the ground state has positive parity. Inspection of Table U shows that the 
parity of the ground state is very sensitive to r± . A small difference in the scale factors r± is more strongly correlated 
with the parity of the ground state than a small difference in the numbers N± of many-body basis states. 

Once the scale factors are determined, we can test how well the right-hand side of Eq. ([T|) can be used to determine 
the parity of the ground state. Our results show that the application of Eq. ([T]) with an energy-independent scale factor 
does not yield reliable predictions. Indeed, Fig. [T] suggests that a linear relation r±(i?ground(i)) = a± + &±£'ground(±) 
should describe the data more accurately. Again, we determine the coefficients a± and b± by fit, and then employ 



8 




FIG. 1: (Color online) Scale factors r± of Eq. ([T]) for a system of m = 9 fermions on = 9 single-particle orbitals with positive 
parity and £2 ~ 9 orbitals with negative parity as a function of the ground-state energies -Eground(±)- 



the right-hand side of Eq. ([T]) with the energy-dependent scale factor to determine the parity of the ground state as 

/ uTtHV- — a^a^ tiTtHV+ — a+cr^ 



sign 



1 + 5_cr_ 



1 + b+a+ 



(15) 



Though this estimate is not correct for each individual member of the ensemble, it yields reasonably reliable predictions 
for the estimated probability p-|_(est) of finding a ground-state with positive parity. Our results for this probability 
are shown in the last column of Table H] 
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TABLE I: Results of numerical simulations. Here, m, £1, and £2 denote the number of fermions and the number of single-particle 
levels with positive and negative parity, respectively. N± is the number of many-body states with the indicated parity, and r± 
denote the scale factors. p+ denotes the probability that the ground state has positive parity, while p+(est) is the probability 
that the estimated ground state has positive parity (based on Eq. ^ with a scale factor that is a polynomial of degree one in 
the energy). 



VI. DILUTE LIMIT 



In canonical random-matrix theory, attention is usually focused on the limit of large matrix dimension. We follow 
suit by considering our model in the "dilute limit" defined by £1^2, m 00 and m/£i^2 0. In practice, we 
compute the leading order of expressions of interest under the strong conditions 1 ^ m ^ £1^2- We show that the 
weight factors appearing in the traces of with k = 1,2 for positive and for negative parity become asymptotically 
equal. That statement holds not only for £1 = £2 but also for £1 ^ £2- 
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Equations ((T0|) and IM]) show that for the positive parity states, all weight factors have the general form 

with ai, q;2, /?2 small positive integers. We evaluate the sums in Eq. (|16p and the corresponding sums defining 
with the help of Stirling's formula, n\ k, expjnlnn — n}. With /i integer, we write TO2 = 2/i, m\ = m ^ 2ii and have 
for the numerator of Eq. (all terms except for (Af+)^^) 

^exp i (^1 - ai)ln(£i - ai) - (m - 2/i - /3i)ln(m - 2/i - 

-(^1 - ai - m + 2/i + ln(^i - ai - m + 2/i + 
+ (^2 - "2) ln(4 - a2) - (2/i - /32) ln(2/i - /32) 



-(£2-a2-2/i + ^2)ln(^2-a2-2/i + /32)| . (17) 

We write the sum as an integral over /i. The integrand assumes its maximum value at 

(0) _ 1 (m - Pi){£2 - 0:2) + f-hjti - ai) . . 

''^ "2 4-«i+£2-«2 ■ ^ ' 

With (5/i = /i — /io, expansion around the maximum yields the negative-definite quadratic form 

2(5/i)2 2((5/i)2 2((5/i)2 



m - 2/i^"^ - /3i ^1 - ai - m + 2/^^°^ + /3i 2/i^^°^ - (32 

2{5y^f 1 (^/.)^ 

^2 - a2 - 2/1^') + /?2 2 r2 



(19) 



Here the last equation defines the width t. For ^ 1, ^2 ^ 1, ^ 1 we have /xp ^ 1. For the dilute limit, we 
neglect terms of higher order, and the resulting integral is Gaussian. We extend the integration from — c» to +00. 
The numerator of expression (jl6p becomes 



'27rr exp |(£i — ai) \n{li — ai) + (^2 — 02) ln(£2 — 012) 
X exp I - (m - 2/i^"' - I3i)\n{m - 2/i^°^ - /3i) 
X exp I - {2^lf - P2) ln(2/iV°^ - /^a) 

X exp I - (€1 - ai - m + 2/i^"^ + /3i) ln(^i - ai - m + 2/i^"^ + (3i) 

X exp ( - {£2 -Oi2- 2/iV"' + ^2) ln(£2 - "2 - 2//°^ + (^2)) ■ (20) 



Using the same approximations to calculate N-^^ we obtain a result of the form (|20p but with ai, a2, /?2 everywhere 

(including the definitions of r and /x^^'') replaced by zero. 

We turn to the negative-parity states. For these states, the word "even" in Eq. ([16]) is replaced by "odd" and iV+ 
by A^_. The calculation is completely analogous except for the replacements (i\ ^ (i\ + \ and /?2 ^ /32 — 1- For the 
maximum of the integrand, that implies 2[i'^ = 2//^'' — 1. As a consequence, the terms 2/i(*'^ -f /?i and 2/i'^°) — (52 
have the same values for states with positive and with negative parity. This in turn implies that the widths r and 
the terms in the exponential in expression (j20p have the same values for states with positive and with negative parity. 
It follows that in our approximation every weight factor for states with positive parity has the same value as the 
corresponding weight factor for states with negative parity. This result is valid beyond the Gaussian approximation 
used in obtaining Eq. (j20p . Indeed, the fundamental form (jl7p depends on [i only through the invariant combinations 
fi + Pi and fi — (32- Modifications can arise only in cases where the limits of integration (which depend on ai,a2, Pi, 
and P2) play a role, i.e., for small values of £1, £2, or m. 
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We have shown that in the dilute hmit and for every reahzation of our random-matrix model, both the first and the 
second moments of H coincide in leading order for states with positive and for states with negative parity. The same is 
true of the matrix dimensions iV+ and N^. Thus for every realization, our Eqs. ([T]) to ([3]) predict equal values for the 
ground-state energies for both parities. How reliable is that prediction? We recall that in the dilute limit, the average 
spectrum of the embedded random two-body ensemble (EG0E(2)) is Gaussian 5]. The proof given in Ref. ^ applies 
likewise to our model. We expect, therefore, that in the dilute limit and to a very high degree of approximation the 
spectrum of any given realization of the ensemble has also Gaussian shape. (For a single realization, the shape of the 
spectrum is defined by taking local averages over a number n N± of neighboring levels.) That expectation rests 
on the plausible assumption that our random-matrix model is ergodic, at least in the dilute limit, and implies that 
for every realization, our Eqs. ([1]) to ([3]) become even better approximations as the matrix dimension increases. We 
conclude that the probabilities for ground states of positive and for negative parity are equal in the dilute limit. That 
conclusion holds with the following proviso. A preference for ground states of, say, positive parity might occur if the 
local spectral fluctuation properties of the two ensembles were locked in such a way that the positive-parity ground 
state fluctuates more often towards smaller energies than does its opposite number. In the next section we exclude 
that possibility. We do so by investigating higher moments of H . 



VII. SPECTRAL FLUCTUATIONS 



Given the coincidence of both the first and second moments of H for states of either parity in the dilute limit, 
we ask: does that coincidence extend to all higher moments so that the local spectral fluctuation properties of both 
ensembles are completely locked? We approach the answer by studying higher moments of H . 

We consider nTrace(7?'^'P±) for k integer and A: > 3. These traces are now shown to have the same structure as 
the first and second moments of H: each trace is a sum of terms each of which is the product of a monomial (or 
polynomial) of order k in the two-body matrix elements (the same for the projectors 7-"+ and V-) and a weight factor 
that does not depend on the random variables but may have a different value for positive and negative parity. 

We proceed as in Section IIVI but are interested only in the general form of the result. The operator is a 
monomial of order k in the matrix elements V^^\V'''^\ X'''^\ X^'^\ Each matrix element carries four indices. Thus, 
in there occur 4fc independent summations over single-particle level indices. Non-vanishing contributions to the 
trace of arise only from Wick-contracted terms. Each pairwise Wick contraction of a creation and an annihilation 
operator in produces a factor of the form ni^, (1 — nic,), n2ct, or (1 — n2a): as the case may be. At the same time, 
two summation indices become equal. After all Wick contractions are done, contains at most 2k independent 
summations over level indices. (That number may be smaller than 2k since two or more of the resulting factors nia, 
(1 — nia), n2a, or (1 — n2a) may carry the same index.) By using the identity — n for the number operator, the 
Wick-contracted H'' can be written in such a way that the summation indices on all such factors are different. For 
k = 2, that was done in Eq. p^ . We consider a single term resulting from this procedure and denote by fci, fc2, fca, ^4 
the powers of the four types of factors (in the same sequence as listed above) in that term. The maximum power with 
which all factors jointly can appear, is 2k so that ki + k2 + k^ + k^ < 2k. Clearly we must also have ki + k2 < (1 and 
^3 + ^4 < ^2 • We conclude that a general term in the Wick-contracted form of H'' , characterized by the four integers 
^1,^2,^3, ^4 as constrained above, has the form 

E E E E 

ai,a2,...,afci /3i,/32,--.,/3fc2 7i i72 ,7fc3 <5i ,^2 ,5fc4 
/■ ki k2 ks k4 

\ n "i"- n ~ ^^1^- ) n "27* n ^ "25„ ) 

^ r=l s=l t=l u=l 

/qi, ...,Qfcj^;/3i,...,/3fc2;7i,...,7t3;(5i, ...... .fc^ ■ (21) 

The sums in this expression are jointly constrained by the condition that no two summation indices are equal. The 
form of the function / depends upon the value of k. / is a monomial of order k in the matrix elements V^, , , X^ . 
These carry the summation indices. The Wick-contraction of H*' yields a sum of terms of the form (I2ip . For the 
calculation of nTrace[i/'''P±], we observe that the expression 



, fci k2 fea ki 

n±(fci, fc2, A:3,/c4) = nTrace<^ J| J|(l - Jt-i/jJ J]^ ^^274 H*^-^ ~ "25„)^± 



(22) 
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does not depend on the values of the indices ai, . . . ,Sk^. Therefore, the normahzed traces of the projections of the 
expression (|2ip are given by 

/Ql,...,afei;/3i,...,/3fc2;7i,...,7fc3;(5i,...,...fc^ 

0!i,...,Qfej /3l,...,,3fc2 7i.--->7f:3 i5i,...,(5fc4 

xn±(fci,fc2,A:3,A:4) . (23) 

Expression ((23)) shows that the results derived in Sect ion HVl for nTrace[if'P±] and for nTra.ce[H^'P±] hold for arbitrary 
powers k of H: Each trace nTrace[iJ'=P±] is a sum of terms; every term in the sum is the product of two factors. The 
first factor contains the random variables and is the same for the states with positive and with negative parity. The 
second factor, a weight factor, may depend on parity. We have, thus, shown that the Hamiltonians for states with 
positive and with negative parity are very highly correlated. 

This is a remarkable result in its own right. Indeed, with increasing values of and £2 the matrix dimensions 
A^-i- and iV_ grow approximately like ((£1 + ^2)/^)™ while the number of two-body matrix elements only grows like 
{£1 +£2)'^- Thus, for {£1 +£2) > Tn? , the matrix dimensions become asymptotically very much larger than the number 
of independent matrix elements. Still, in the sense of Eq. (|23p . the two Hamiltonians remain totally correlated. 

We turn to the weight factors appearing in Eq. ([23]) and show that these are also asymptotically equal. Our statement 
applies up to a maximum value of k which we determine approximately. The weight factors n±(A:i, A:2, fcs, ^4) are 
explicitly given by 

n+(fci, /C2, /Ca, fc4) = — — (5mi+m,.m^m2,ovon 1 ^ 7^ ) f ^ 7'' ) (24) 

N+ ^ ' \mi - kij \m2 - k^J 



mi .ni2 



and 



n_(A:i, ^2, fcs, /C4) = — — <^mi+m2,m<5m2,odd 1 ^ ,^ M ^ 7^ ) ■ (25) 

V"^i - klj \m2 ~ k^J 



nil ,m2 



In the summations over mi,TO2, we obviously must have mi > ki and 7712 > k^. Since mi and TO2 are both bounded 
by m, that condition in fact limits ki and k^. It is obvious that for large values of fc, the two weight factors cannot 
always be equal. Consider, for instance, the case fci = 0, k^ — m. Then we have mi ~ and 777,2 = w. That implies 
n+(fci, k2, ks, k4) = 0, n_(fci, fe, fca, ki) ^ if m is odd and n_(fci, k2, fcs, ki) = 0, n+(fci, k2, fcs, ki) ^ if 777 is even. 
To avoid such cases, we must have k < m. Even then 11+ and n_ may differ. This happens when the bounds on the 
summation indices in Eqs. ([M]) and pS)) are relevant. We avoid these cases by choosing fc <C to. We recall that the 
asymptotic regime is characterized by the relations 1 ^ 777, <C ^1 , ^2 ■ We thus require that 777 is sufficiently large to 
accommodate the relation k <^ m and yet allows k to assume values large compared to one. With these assumptions, 
the arguments used above for k — 1,2 show that 11+ = n_. 

We have shown that in the asymptotic regime and for all k with k < kg, the moments Trace {H'^'P±) pairwise have 
the same values for states with positive and with negative parity. Here fco obeys 1 ^ fcg <SC 777. That conclusion does 
not depend on assuming any symmetry such as ii = £2- We have also shown that for k ko, the moments differ. 
As k increases, the bounds on the summations over products of binomial factors become ever more important. As 
a consequence, the differences between moments for states with positive and with negative parity increase with k. 
That statement is relevant for the local spectral fluctuation properties of both ensembles. Indeed, it is known [8] that 
such fluctuation properties depend on the very highest moments of H: In the limit of infinite matrix dimension, there 
exists a clear separation between the overall shape of the spectrum (defined by averaging over an energy interval large 
compared to the average level spacing d), and the local spectral fluctuations (defined on a scale of order d). Since 
the moments of H for states of positive and negative parity differ for fc 3> fco , we conclude that the local fluctuation 
properties of both ensembles are uncorrelated in the dilute limit, even though the moments of H for both parities 
coincide up to fc « fco . This excludes the possibility mentioned in Section IVII that the local spectral fluctuation 
properties of the two ensembles are locked in such a way that the positive-parity ground state fluctuates more often 
towards smaller energies than does its opposite number or vice versa and completes the proof that in the dilute limit, 
ground states of either parity carry equal probabilities. 



VIII. SUMMARY AND DISCUSSION 



We have shown that in the dilute limit, ground states of either parity carry equal probabilities. That conclusion is 
based on the following facts, (i) The spectra are asymptotically Gaussian, and Eqs. ([T]) to ([3]) become asymptotically 
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strictly valid, (ii) The first and second moments of H and the dimensions of the Hamiltonian matrices become 
asymptotically equal for either parity so that Eqs. H]) to ([3]) predict equal probabilities for either parity, (iii) The 
local spectral fluctuation properties of the two spectra are asymptotically uncorrelated because very high moments 
have different values. That fact excludes a locking of these fluctuations. 

Deviations from equal ground-state probabilities are, thus, finite-size effects. For values of the parameters to, £1,^2 
that are sufficiently small for numerical simulations, we have indeed found such deviations. They occur whenever 
the dimensions and iV_ differ. Conversely, for Nj^. = 7V_ we have not found significant deviations from equal 
probabilities. The small fiuctuations found for r± in the fits to the data show that Eqs. ((T|) and ^ are approximately 
valid: They do predict correctly which parity has the higher probability to furnish the ground state. The values of 
the predicted probabilities are semi-quantitatively correct. 

Calculations using the two-body random ensemble (TBRE) reported in Rcf. [10] displayed correlations between 
spectra carrying different quantum numbers. One may argue that these results contradict our present findings. This is 
not the case: Calculations using the TBRE are neccessarily restricted to small matrix dimensions while our argument 
for independence of spectral fluctuation properties of states of positive and negative parity applies only in the dilute 
limit, i.e. for infinite matrix dimension. 

Our results may have interesting implications for the statistical theory of nuclear reactions. There an open question 
is this: are ^-matrix elements carrying different quantum numbers like total spin uncorrelated? That assumption is 
always used in the theory and is consistent with the observed symmetry of compound-nucleus cross sections about 90 
degrees in the center-of-mass system. Still, the assumption is not obviously valid for a realistic random-matrix model of 
nuclear reactions. Normally the statistical theory of nuclear reactions uses the Gaussian orthogonal ensemble (GOE). 
It would be more realistic to use instead the TBRE. The TBRE differs from the GOE in that it employs a shell-model 
in which the two-body matrix elements are the random variables. (For a review of the TBRE, see Ref. [ll|.) But then 
it is the same set of random variables that govern scattering matrix elements carrying different quantum numbers; 
just as in the model considered above the same random two-body matrix elements govern the Hamiltonians for states 
of different parity. To approach the question, we observe that for orthogonally invariant ensembles, universality holds 
also for elements of the scattering matrix carrying identical quantum numbers [l^ . That statement implies that 
correlations between such elements depend only on local spectral fluctuation properties. This conclusion is supported 
by the explicit calculation in Rcf. [13] of the correlation function of a pair of S'-matrix elements: Aside from the 
strength of the coupling to the open channels, the correlation depends solely on the value of the local mean level 
density. If we assume that these statements carry over to the TBRE, and if we further assume that in the TBRE just 
as in the model studied above, the local spectral fluctuation properties of spectra carrying different quantum numbers 
are uncorrelated in the limit of large matrix dimension, we are led to the conclusion that S'-matrix elements carrying 
different quantum numbers are, likewise, uncorrelated. The limit of a large matrix dimension is appropriate because 
the resonances relevant in the statistical theory correspond to states above the ground state. 
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